Level spacing distribution of a Lorentzian matrix at the spectrum edge 
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Effective Hamiitonians can expiain in a much simpier way the physics behind a scattering process. 
Chaotic scattering is directiy reiated to Lorentzian Hamiitonians which, because of their properties, 
can be reduced to a 2 x 2 matrix problem in the case of two modes scattering. In this framework, 
we provide the distribution of level spacing of its eigenvalues and show that this special kind of 
distribution has no mean level spacing (divergent) and is characterized by a geometrical decay law. 
We discuss the relation of this distribution to the averaged level spacing at the edge of the spectrum 
of TV x N Lorentzian matrices. 

PACS numbers: 05.45.-a, 72.10.-d, 73.23.-b 



(N 



i 

C 

O 

o 



>: 

Ov 

in 



In random matrix theory, the level spacing, defined 
as the distance between two adjacent eigenvalues s — 
\E a — E a ^i\, plays an important role in the study of 
chaotic cavities. Indeed, the distribution of this vari- 
able shows whether the system is integrablc (Poisso- 
nian statistics for the level distribution V{s)) or non- 
integrable (strong level repulsion at small s)[l]. More- 
over, the mean level spacing is an important scaling pa- 
rameter for many physical observalbcs. One of the most 
astonishing success of random matrix theory is the expla- 
nation of the distribution of the eigen-energies in nuclear 
spectra. Wigner obtained this distribution from the level 
spacing of a 2 x 2 Gaussian matrix. Remarkably, it has 
been noticed that even though the resulting distribution 
is not exact for the N x N case [1][2][3], it remains an ex- 
cellent approximation and still referred to as the Wigner 
surmise distribution. In this framework, we will be in- 
terested in the distribution of a 2 x 2 Lorentzian matrix. 
The choice of this case is very subtle and goes beyond 
the simple mathematical curiosity. The first argument 
is related to Dyson's circular ensembles (CE). Indeed, in 
order to obtain a scattering matrix with an "equal a pri- 
ori probability ansatz", we define an ensemble of unitary 
matrices 5* uniformly distributed according to the Haar 
measure P(S)dS = -gSS. C is a normalization constant. 
The scattering matrix S relating the amplitudes of the 
incoming wave functions to the outgoing ones is related 
to the scattering Hamiltonian as follows [ ] : 
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implies uniformly distributed scattering matrices 
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P stands for the different symmetry cases (/3 = 1 in pres- 
ence of time reversal symmetry TRS, /3 — 2 without TRS, 
and j3 = 4 broken spin-rotational symmetry) and N is 
the size of the Hamiltonian. The center £ and the width 
A are directly related to the self energy of the channels 
E = £ — iX. It is worth stopping at this interesting result 
and precise that is exact and valid for all the Hamilto- 
nian sizes N. Nevertheless, for very large systems, it has 
been shown [5] that a Gaussian matrix exhibiting a same 
mean level spacing becomes equivalent to the Lorentzian 
distribution. Actually, this result is more general and 
stands for each distribution P(H) oc exp {— NtrV(H)} 
with a confining potential V generating a smooth mean 
level spacing A ~ 0(N^ 1 ) when /V^oo [6] [7]. It comes 
out from this remark that the distribution of the local 
level spacing of a N x N Lorentzian matrix (N being 
large) is given by the Wigner surmise when the spacing 
is taken in the bulk of the spectrum. The 2x2 Lorentzian 
matrix do not have a dense spectrum and therefore do not 
fulfill this condition. Now, we come back to the second 
argument concerning the choice of a 2 x 2 matrix. The 
scattering of two independent and equivalent conducting 
modes is widely considered. In this case, the size of the 
scattering matrix is 2x2. Nevertheless, the Hamiltonian 
is an N x N matrix where the size N is generally taken 
very large. Since there is only two modes, we can rewrite 
the scattering matrix as follows [8] [9] : 

(3) 



S = 1 - 27riVV t - 



The N x N matrix % is the Hamiltonian of the chaotic 

cavity and £ is the self energy of the channels (the con- _ %(E) — £ 

ducting modes). W is a coupling matrix between the 
leads (channels) and the Hamiltonian of the closed cavity. 
Thus, for a finite size Hamiltonian describing a chaotic 
cavity, It can be shown [5] that considering a Lorentzian 

H{E) 

distribution for %, with appropriate center and width [8], taining only the two conducting modes. £ is a 2 x 2 



where "H(£0 is an effective Hamiltonian, energy depen- 
dent, which has for size 2x2. It is defined as follows 
e-u t - W is a 2 x 2 sub-matrix of W con- 



2 



submatrix of £ containing only the relevant terms, r is 
N x 2 matrix with rj- = Sij. This operation is exact 
and not an approximation (up to a phase depending on 
the choice of the scattering potential limits). It allows to 
reduce consequently the high degrees of freedom of the 
original cavity. The amazing thing is that the resulting 
effective Hamiltonian % is still Lorentzian. This is due 
to the following two properties [5] [10]: 

• If A is Lorentzian with center and width (£ , A), 
is also Lorentzian with ( g2 ^ A2 , £2 ^ A2 ) . 

• If A is Lorentzian, any submatrix of A is Lorentzian 
with the same width and center. 

To summarize, we can say that any N x N Lorentzian 
Hamiltonian describing the cavity, boils down to a 2 x 2 
effective Lorentzian matrix problem. This, suggests that 
the transport coefficients (transmission, shot noise . . . ) 
are closely related to the statistics of the two levels of the 
effective Hamiltonian H(E). In what follows, we give the 
distribution of the level spacing of this 2x2 Lorentzian 
matrix. 



Level spacing distribution 

The starting point to express the level spacing of a 2 x 2 
Lorentzian matrix % is the distribution of the eigenval- 
ues (we call them x and y). This is obtained from the 
distribution Eq. (2) by the usual variable change of the 
eigenvalues [5] [12]. We start with the case of a centered 
{£ = 0) and standard (A = 1) distribution: 

\x-y\ p 



P{x,y) 
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(4) 



27T (1 + Z 2 )W2+1)(1 + 2,2)03/2+1) 

I a; — y\P is the Jacobian of the variable substitution. It 
acts as a repulsion term avoiding the eigenvalues to be 
the same [1][12]. The distribution of the level spacing s 
is therefore expressed as follows: 



V(s) 



S(s -\x-y\ 
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dxdy 



[(l + .T 2 )(l + y 2 )]W/ 2 + 1 ) 2tt 



After a first integration over the delta function, one ob- 
tains: 



(5) 



= — lp(s)+p(-s)} 
where the function p(s) reads: 

( ] = r +oc dx 

P[S> J_ 00 [l + x 3 ]W 3 + 1 )[l + (a-x) 2 ]W 3 + 1 

It is interesting to notice that p(s) = (/ * /)(s) 

where * is the convolution product and / is the following 

function: 
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m (6) 



f{x) = 



( 1 + 3-2)03/2+1)' 



(7) 



The convolution product suggests to use the Fourier 
transform ^(notcd also " ) : 



Hp{s)) = f 2 



(8) 



The Fourier transform of / Eq. (7) is expressed using 
the modified Bcsscl function of the second kind, BesselK 
[11]: 



/(w) oc | w |^BesselK[^ti, | w [ 



(9) 



To obtain the function p(s), we need to do the inverse of 
the Fourier transform of Eq. (8). 



p(s) (x J- 1 ( (M^BesselK[^, M]) 2 ) 



The inverse Fourier transform in Eq. (10) is expressed 
using the Gauss hypcrgcometric function 2-F1 [11]: 



p(s) oc 2 F 1 (- 



fi + 2 (3 + 3 
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Now, we are ready to express the probability density 
function of the level spacing for the three ensembles in a 
unique compact form deduced directly from Eqs.(5) and 
(11): 

V(s) = C^ 2 F 1 ( 3 -+P,^,^,-^) (12) 
where the normalization constant Cp is given as follows: 

1 r(I + /j) 



a 
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(13) 



The behavior of this probability density at small level 
spacing (s — > 0) is given by: 



(14) 



this is the signature of the level repulsion of each of the 
three ensembles (/3 = 1, 2 or 4) corresponding to different 
type of symmetry. This feature is comparable to what is 
found for the Wigner surmise[l][12]. The tail of the dis- 
tribution is more interesting since it is different from all 
the distributions known so far in a sense it has a geomet- 
rical fall off: 



(15) 



this tail is /3-independent and the exponent of the fall 
off law implies the important property of a no mean for 
the density function. It means that the mean level spac- 
ing defined usually as s = J sV(s)ds does not exist since 
the integral diverges. This signifies that despite the dis- 
tribution is narrow around it's center (for small A), the 
eigenvalues can be typically as far as possible. 



THE DISTRIBUTION OF THE LEVEL SPACING 
AT ARBITRARY WIDTH 

It is obvious that changing the center of the distribu- 
tion does only shift all the energies by the same amount 
and therefore the level spacing distribution remains un- 
changed. This is not the case of the distribution width . 
Indeed, at a different energy, the width of the Lorentzian 
changes and becomes arbitrary A. We can obtain the 
distribution at this energy by changing the variables in 
the integral we started with or just by saying that with 
this new width of the distribution, all the eigenvalues are 
multiplied by a factor A and therefore, the invariant dis- 
tribution is obtained for the renormalized level spacing 
defined as follows: 

*-j 

Within this definition, the distribution of the normalized 
variable S is the same as in Eq. (12). Usually, we prefer 
to normalize the lengths in the problem with the mean 
level spacing s. Here, since the mean level spacing is 
divergent, we can not use it for this task. There is an- 
other quantity which can be used as a scaling length. It 
is also commonly called in literature mean level spacing 
and noted A. This variable is defined as the inverse of 
the density of states at the center of the Hamiltonian 
spectrum: A = p ^e=s) wncrc the density p is defined 

as p{E) = T,E i ={x,y}( 6 ( E - R i))- This scale A is sti11 
defined for the case of a 2 x 2 matrix and is finite. Nev- 
ertheless, it has no anymore the signification of a mean 
level spacing since the spectrum is not dense. The den- 
sity of states p is j3- independent and reads [5] : 



P(E) 



A 



(16) 



ir{E- E) 2 + A 2 
thus, the scale A at the center of the distribution reads: 

(17) 



A = ^A 
2 



This parameter is equivalent (up to a factor §) to what 
was stated before about choosing the width A as a simpler 
and natural scaling parameter. Moreover, by choosing A, 
the factors in Eq. (12) remain unchanged. Hereafter, the 
variable s will stand for the level spacing renormalized 
by the width of the distribution A. Within this defini- 
tion, it is worth keeping in mind the important feature 
of the level spacing distribution at small s (resistance to 
crossing) and at large s (no mean level spacing) which 
are true for the three ensembles ((3 = 1, 2 and 4). 




if s < 1 
if s > 1 



(18) 



We stress that this geometrical fall-off exhibits a high de- 
gree of fluctuation of the level spacing. It is completely 
different from the usual decay we face in the Wigner sur- 
mise or the semi-poissonian ensembles[13][14]. 




FIG. 1. The distribution of the renormalized level spacing 
of a 2 x 2 Lorentzian orthogonal matrix (/? = 1). The red 
line is the analytical result given in Eq. (19). The agreement 
between theory and numerical simulation is excellent. 



Orthogonal case /3 = 1 

It may be interesting to simplify formula Eq. (12) for 
each case of the three symmetries (/? = 1,2 and 4). We 
start here with the case of time reversal symmetry j3 = 1. 
It can be shown that the expression of the level spacing 
distribution boils down to the following form using the 
complete elliptic integrals of the first and the second kind 
noted respectively K and E [111: 



V{s) = 



(-4 + s 2 )E(-^) + (4 + s 2 )K(-^) 



(19) 



s(4 + s 2 ) 2 

The numerical simulation consisting on sampling 2x2 
Hamiltonians with a Lorentzian distribution and doing 
the statistics of the level spacing shows an excellent agree- 
ment with formula Eq. (19) as can be shown in Fig. 1. 
The expansion at small and large level spacing can be 
obtained straightforwardly: 



V(s) 




4r, ifs<l 



if s > 1 



UNITARY CASE /3 = 2 

The case of the unitary symmetry (/? = 2) can be writ- 
ten in a much simpler formula. Just by taking into ac- 
count the form of the Gauss hypcrgcometric function for 
even parameter (/?) , one finds the following result: 



V{s) = 



4 s 2 (20 + s 2 ) 
n (4 + s 2 ) 3 



(20) 



The behavior at small and large level spacing is in agree- 
ment with formula Eq. (18): 



T{s) 




if s < 1 
if s > 1 
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FIG. 2. The distribution of the renormalized level spacing of 
a 2 x 2 Lorentzian Unitary matrix (/3 = 2). The red line is the 
analytical result given in Eq. (20). The agreement between 
theory and simulation is excellent. 



P(s) 




FIG. 3. The distribution of the renormalized level spacing 
of a 4 x 4 Lorentzian matrix with the symplectic symmetry 
(/3 = 4). (notice the factor 2 due to the degeneracy of the 
eigenvalues). The red line is the analytical result given by 
Eq. (21). The agreement between theory and simulation is 
excellent. 

This distribution law is tested numerically (see Fig. 2). 
Again, we notice that this distribution has no mean level 
spacing. 



/? 

FIG. 4. The exponent of the geometrical fall-off l/x a of an 
N x N Lorentzian matrix as a function of the number of levels 
n = N e taken to define the averaged mean level spacing. The 
levels are taken at the edge of the spectrum starting from the 
largest eigenvalue. 
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FIG. 5. Distribution of the averaged level spacing at the edge 
of the spectrum for N e = 8 and j3 = 1 (Blue curve). The 
size of the matrix is N = 40. The red curve represents the 
theoretical result for 2x2 orthogonal Lorentzian ensembles. 
The figure shows a very good agreement in the description of 
the tail which suggests a very high degree of fluctuation and 
the absence of the mean for this distribution. 



Symplectic ensemble /? = 4 

Again, the distribution of the level spacing for the sym- 
plectic ensemble is much simpler and can be written in a 
fractional form: 



Now, it becomes clear that the tail of the level spacing 
distribution is j3- independent. 

The numerical simulation showing an excellent agree- 
ment with the distribution Eq. (21) is given in Fig. 3. 



N x N Lorentzian Hamiltonian 



4s 4 (336 + 24s 2 + s 4 ) ^ * s WG ^ known ["'] * na * a iar g e matrix taken from the 

V(s) — — (21) Lorentzian ensemble is equivalent to a Gaussian matrix 



(4 + s 2 ) E 



sharing the same mean level spacing A at the bulk of 



The expansion at small and large level spacing is straight- the spectrum. This result was obtained by comparing 
forward: the cluster functions [5] and one can understand this as 

a special case of the general idea of the bulk spectrum 
universality [6]. Therefore, it is easy to test that the 
level spacing taken from the bulk spectrum is well fitted 
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FIG. 6. log(V(s)) as a function of log(s). The linear decay 
for large values of s is characteristic of a geometrical fall-off 
of the distribution V(s). Red curve is the linear fit of the 
numerical data. It has for equation: y — —1.95a; + 0.30. The 
absolute value of the slope a = 1.95 represents the exponent 
of the geometrical fall-off. 



by the Wigner surmise. The situation is not the same 
at the edge (tail) of the spectrum: the fluctuations are 
higher and the mean level spacing diverges. Fig. 5 shows 
the distribution of the averaged level spacing of an Or- 
thogonal Lorentzian matrix at the edge of the spectrum, 
defined as : 



1 



i=l 



E{ — Ei 



where the sum is taken over some iV e levels at the edge of 
the spectrum (the levels Ei are ordered such that E\ rep- 
resents the largest eigenvalue). It shows clearly that the 
tail falls off with a geometrical law comparable to that of 
the 2x2 matrix studied in the previous sections. This 
can be more visible in the plot of log('P) as a function of 
log(s) shown in Fig. 6. The decay law seems quite well 
fitted by 4r with a a parameter depending on the num- 
ber of levels considered in the edge. As shown in Fig. 5, 
the law of the averaged level spacing has no mean since 
mostly we have a < 2. It is now important to define 
where the edge of the spectrum roughly starts. This can 
be seen as the levels where the distribution of the spacing 
between two consecutive levels stops to be well fitted by 
a Wigner surmise. In Fig. 4, we give the exponent a 
for different numbers N e of levels taken from the edge of 
the spectrum. We see that the exponent a approaches 
the value a — 2 when we increase N e . We stress that 
increasing N e needs larger values of N in order to count 
only the eigenvalues at the edge. 

This behavior at the edge of the spectrum suggests a new 
physics different from that of the bulk. This was already 
noticed in [8] [9] where the statistics of thermopower and 
the delay time of a chaotic cavity were found to be differ- 
ent in the two situations corresponding to a Fermi energy 



lying in the bulk of the spectrum or in its edge. Moreover, 
some features and distributions can be directly obtained 
from the 2x2 Lorentzian Hamiltonian instead of the 
original N x N matrix. Indeed, the distribution of the 
Seebeck coefficient [8] and the Wigner's time [9] at the 
edge of the spectrum were found directly by considering 
the 2x2 Lorentzian Hamiltonian. 



Conclusion: 

We gave the exact form of the level spacing distribu- 
tion of a 2 x 2 Lorentzian matrix. This kind of matrices 
appear in chaotic scattering problems as an effective ma- 
trix replacing large matrices describing chaotic cavities. 
The tail of the distribution is slowly decaying and leads 
to the absence of the mean level spacing. 
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